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Abstract 
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ters it is possible to write the new series solutions in terms of elementary functions. We 
regain well known physically reasonable models. A physical analysis indicates that the 
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field on the gravitational interaction is highlighted. 
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1 Introduction 



Solutions of the Einstein-Maxwell system of equations for static spherically symmetric 
interior spacetimes are necessary to describe charged compact objects in relativistic 
astrophysics where the gravitational field is strong as in the case of neutron stars. In 
the presence of an electromagnetic field, the gravitational collapse of a spherically sym- 
metric distribution of matter to a point singularity may be avoided: the gravitational 
attraction is counterbalanced by the Coulombian repulsive force in addition to the 
pressure gradient. The recent analyses of Ivanov [1] and Sharma et al [2] show that 
the presence of the electromagnetic field affects the values of redships, luminosities and 
maximum mass of a compact relativistic object. Patel and Koppar [3], Tikekar and 
Singh [4], Mukherjee [5] and Gupta and Kumar [6] demonstrated that it is possible to 
model charged neutron stars with high densities with acceptable bounds for the surface 
redshift, luminosity and total mass. Thomas et al [7], Tikekar and Thomas [8] and Paul 
and Tikekar [9] demonstrated that charged relativistic solutions may be applied to core- 
envelope models. The role of the electromagnetic field in describing the gravitational 
behaviour of stars composed of quark matter has been recently highlighted by Mak 
and Harko [10] and Komathiraj and Maharaj [11]. Therefore, the Einstein-Maxwell 
system for a charged star has attracted considerable attention in various physical in- 
vestigations. We note that at present there is no unified theory of electromagnetism 
and gravitation which is generally accepted. In this paper we have used the approach 
of coupling the electromagnetic field tensor to the matter tensor in Einstein's equations 
such that Maxwell's equations are satisfied. We believe that the qualitative features 
generated in this charged model should yield results which are physically reasonable. 

Our objective is to generate a new class of solutions to the Einstein-Maxwell system 
that satisfies the physical criteria: the gravitational potentials, electric field intensity 
and matter variables must be finite and continuous throughout the stellar interior. The 
speed of the sound must be less than the speed of the light, and ideally the solution 
should be stable with respect to radial perturbations. A barotropic equation of state, 
linking the isotropic pressure to the energy density, is often assumed to constrain the 
matter distribution. In addition to these conditions the interior solution must match 
smoothly at the boundary of the stellar object with the Reissner-Nordstrom exterior 
spacetime. In recent years researchers have attempted to introduce a systematic ap- 
proach to finding solutions to the field equations. Maharaj and Leach [12] generalised 
the Tikekar superdense star, Thirukkanesh and Maharaj [13] generalised the Durgapal 
and Bannerji neutron star and Maharaj and Thirukkanesh [14] generalised the John 
and Maharaj [15] model. These new classes of models were obtained by reducing the 
condition of pressure isotropy to a recurrence relation with real and rational coefficients 
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which could be solved by mathematical induction, leading to new mathematical and 
physical insights in the Einstein-Maxwell field equations. We attempt to perform a 
similar analysis here to the coupled Einstein-Maxwell equations for a general form of 
the gravitational potentials with charged matter. We find that the generalised con- 
dition of pressure isotropy leads to a new recurrence relation which can be solved in 
general. 

In this paper, we seek new exact solutions to the Einstein-Maxwell field equations, 
using the systematic series analysis, which may be used to describe the interior rel- 
ativistic sphere. Our objective is to obtain a general class of exact solutions which 
contains previously known models as particular cases. This approach produces a num- 
ber of difference equations, which we demonstrate can be solved explicitly from first 
principles. We first express the Einstein-Maxwell system of equations for static spheri- 
cally symmetric line element as an equivalent system using the Durgapal and Bannerji 
[16] transformation in Section 2. In Section 3, we choose particular forms for one of 
the gravitational potentials and the electric field intensity, which reduce the condition 
of pressure isotropy to a linear second order equation in the remaining gravitational 
potential. We integrate this generalised condition of isotropy equation using Frobenius 
method in Section 4. In general the solution will be given in terms of special functions. 
However elementary functions are regainable, and in Section 5, we find two category of 
solutions in terms of elementary functions by placing certain restriction on the param- 
eters. In Section 6, we regain known charged Einstein-Maxwell models and uncharged 
Einstein models from our general class of models. In Section 7, we discuss the physical 
features, plot the matter variables and show that our models are physically reasonable. 
Finally we summarise the results found in this paper in Section 8. 

2 The Field Equations 

The gravitational field should be static and spherically symmetric for describing the 
internal structure of a dense compact relativistic sphere which is charged. For describ- 
ing such a configuration, we utilise coordinates (x a ) = (t, r, 9, 0), such that the generic 
form of the line element is given by 

ds * = _ e 2u(r) dt 2 + e 2X(r) dr 2 + ^2 + ^2 M(j) 2y (1) 
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The Einstein field equations can be written in the form 

l[r(l-e^)]' = p (2a) 

4( i - e " 2A )+v :e " 2A = p (2b) 

e - 2X (u" + u' 2 + --u'X'--) = p (2c) 
\ r r J 

for neutral perfect fluids. The energy density p and the pressure p are measured relative 
to the comoving fluid 4-velocity u a = c~ v 5q and primes denote differentiation with 
respect to the radial coordinate r. In the system (I2aj) - (l2"cl) . we are using units where 
the coupling constant = 1 and the speed of light c = 1. This system of equations 
determines the behaviour of the gravitational field for a neutral perfect fluid source. 
A different but equivalent form of the field equations can be found if we introduce the 
transformation 

x = Cr\ Z(x) = e" 2A(r) , A 2 y 2 (x) = e 2iy(r) (3) 
so that the line element ([1]) becomes 

ds 2 = -A 2 y 2 dt 2 + -^-dx 2 + ^-{d9 2 + sin 2 Odcf 2 ). 
ACxZ C 

The parameters A and C are arbitrary constants in Under the transformation (jHJ), 
the system (T2a l) -f l2cl) has the equivalent form 

i^- 2 Z = L (4a) 

izi + — = ^ (4b) 
y x C 

AZx 2 y + 2Zx 2 y + (Zx-Z + l)y = (4c) 

where dots represent differentiation with respect to x. 

In the presence of an electromagnetic field the system (j4"a|) - (j4"c|) has to be replaced 
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by the Einstein-Maxwell system of equations. We generate the system 

- M 

( E 2 x\ 
AZx 2 y + 2Zx 2 y + ( Zx - Z + 1 - — jy = (5c) 

a 2 AZ / ■ \2 

- = + (5d) 

where i£ is the electric field intensity and a is the charge density. This system of 
equations governs the behaviour of the gravitational field for a charged perfect fluid 
source. When E = the Einstein-Maxwell equations (I5a| - (l5d| ) reduce to the uncharged 
Einstein equations fj4aj) - fHc|) . 



3 Choosing Z and E 

We seek solutions to the Einstein-Maxwell field equations (15ap - (15dp by making explicit 
choices for the gravitational potential Z and the electric field intensity E on physical 
grounds. The system fl5aj) - f!5dl) comprises four equations in six unknowns Z,y, p,p, E 
and a. Equation ( l5cl) . called the generalised condition of pressure isotropy, is the master 
equation in the integration process. In this treatment we specify the gravitational 
potential Z and electric field intensity E, so that it is possible to integrate fl5c|) . The 
explicit solution of the Einstein-Maxwell system (I5aj) - (15dj ) then follows. We make the 
particular choice 

Z = \±^ (6) 
1 + bx w 

where a and b are real constants. The function Z is regular at the centre and well 
behaved in the stellar interior for a wide range of values of a and b. It is important to 
note that the choice ([6]) for Z is is physically reasonable. This form for the potential 
Z contains special cases which correspond to neutron star models, eg. when a = — | 
and b = 1 we regain the uncharged dense neutron star of Durgapal and Bannerji [16]. 
When a is arbitrary and 6=1 then Thirukkanesh and Maharaj [13] and Maharaj 
and Komathiraj [17] found charged solutions to the Einstein-Maxwell system. These 
solutions can be used to model a charged relativistic sphere with desirable physical 
properties. Consequently the general form contains known physically acceptable 
uncharged and charged relativistic stars for particular values of a and b. We seek 
to study the Einstein-Maxwell system with the choice §6§ in an attempt to find new 
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solutions, and to show explicitly that cases found previously can be placed into our 
general class of models. 

Upon substituting fl§D in equation (I5c|) we obtain 



4(1 + ax)(l + bx)y + 2 (a - b)y + 



b(b - a) 



E 2 (l + bx) 2 ' ] 



y = o. (7) 



The differential equation (J7J) is difficult to solve in the above form; we first introduce a 
transformation to obtain a more convenient form. We let 

l + bx = X, Y(X)=y(x), b^O. (8) 

With the help of ©, © becomes 



« [o X-(«-6)]^ + 2(a-6)g + 



(b- a) 



E 2 X 2 



Y = (9) 



C(X-l) 

in terms of the new dependent and independent variables Y and X respectively. The 
differential equation (Q may be integrated once the electric field E is given. A variety 
of choices for E is possible but only a few have the desirable features in the stellar 
interior. Note that the particular choice 

2 aC(X - 1) aCbx 

E = — = (TTteF { ] 

where a is a constant has the advantage of simplifying (Q. The electric field given 
in (fIU|) vanishes at the centre of the star, and remains continuous and bounded for 
all interior points in the star. When 6 = 1 then E 2 reduces to the expression in the 
treatment of Maharaj and Komathiraj [17]. Thus the choice for E in ffTO]) is physically 
reasonable in the study of the gravitational behaviour of charged stars. With the help 
of ffTU]) we find that ([HD becomes 

d 2 Y dY 
AX [aX - (a - b)} — + 2(a - &)— + [(6 - a) - a] Y = 0. (11) 

uA fflA 

The differential equation (TIT]) becomes 

d 2 Y dY 
AX [aX -(«-&)] — + 2(o - b)^ + (b-a)Y = (12) 

when a = and there is no charge. 



4 Solutions 

We need to integrate the master equation (fTTj) to solve the Einstein- Maxwell system 
(j5aj) - ([5dj) . Two categories of solution are possible when a = b and a ^ b. 
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4.1 The Case a = b 



When a = b equation ffTTj) becomes 



which is an Euler-Cauchy equation. The solution of (TT3"j) becomes 

Cl (l + ax) (1+ v /I +^)/ 2 + c 2 (l + ax) (1 -v /I +^)/ 2 if a > 0, 



+C2 cos 



' i sin ^ln(l + ax; 
^fln(l + ax; 



if a < 0, 



(13) 



(14) 



where c\ and c 2 are constants. From fl5a|) and ([6]) we observe that p — We do 

not pursue this case to avoid negative energy densities. It is interesting to observe that 
when a = b = then it is possible to generate an exact Einstein-Maxwell solution to 
(|5aj) -( l5dl) . for a different choice of E 2 , which contains the Einstein universe as pointed 
out by Komathiraj and Maharaj [18]. 



4.2 The Case a ^ b 

Observe that it is not possible to express the general solution of the master equation 
(flTl) in terms of conventional elementary functions for all values of a,b (a ^ b) and 
a. In general the solution can be written in terms of special functions. It is necessary 
to express the solution in a simple form so that it is possible to conduct a detailed 
physical analysis. Hence in this section we attempt to obtain a general solution to the 
differential equation ( fTTj) in series form. In a subsequent section we show that it is 
possible to find particular solutions in terms of algebraic functions and polynomials. 

We can utilise the method of Frobenius about X = 0, since this is a regular singular 
point of the differential equation (ITT1) . We write the solution of the differential equation 
(fTT!) in the series form 

oo 

Y = Y,c n X n+r , C ^0 (15) 

n=0 

where c n are the coefficients of the series and r is a constant. For an acceptable solution 
we need to find the coefficients c n as well as the parameter r. On substituting (fT5l) in 
the differential equation fTTTl) we have 

CO 

2(a - b)c r[-2(r - 1) + \}X r - 1 + l 2 ( a ~ b)c n+1 (n + r + l)[-2(n + r) + 1] 

n=l 

+c n [Aa(n + r)(n + r -1) - (a-b + a)}} X n+r = 0. (16) 
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For consistency the coefficients of the various powers of X must vanish in (fl6l) . Equat- 
ing the coefficient of X r ~ x in ffTBl) to zero, we find 

(a - b)c r[2(r - 1) - 1] = 

which is the indicial equation. Since Cq ^ and a ^ b, we must have r = or r = |. 
Equating the coefficient of X n+r in f)16p to zero we obtain 



c n +i 



Aa(n + r)(n + r — 1) — [a — b + a] 



c n , n > 



(17) 



2(a - 6)(n + 1 + r)[2(n + r) - 1] 

The result (iTTj) is the basic difference equation which determines the nature of the 
solution. 

We can establish a general structure for all the coefficients by considering the leading 
terms. We note that the coefficients C\, C2, C3, ... can all be written in terms of the leading 
coefficient c , and this leads to the expression 



n£a(p + r)(p + r — 1) — (a — b + a) 
~/ 7T7 — ~ — ; — \r7T7 — ; — \ — 7~n c o 



p=0 



2(a-b)(p+ l + r)[2(p + r) - 1] 



(18) 



where the symbol Yl denotes multiplication. It is also possible to establish the result 
( fl8l) rigorously by using the principle of mathematical induction. We can now generate 
two linearly independent solutions from ffl5|) and ffl8l) . For the parameter value r = 
we obtain the first solution 



Vi 



Co 



1 + £ n Aap(j) ~ ~ ~ ~ ~ h + a ^ x n+i 



n=0 p=0 



2(o-6)(p+l)(2p- 1) 



(19) 



c 



n=0 p=0 



2(a-6)(p+l)(2p- I] 



For the parameter value r = | we obtain the second solution 



Y 



CqX* 



n=0 p=Q 



a(2p + 3)(2p+l)-(a-b + a) xn 



(+1 



2/2 



c (l + bx) 



(a-b)(2p + 5)(2p + 2) 

1 _u rr q ( 2 p + 3 )( 2 p + jj - (a-fe + q) +1 

i + Z^ll ( a _ 6 )(2p + 5 )(2p + 2) 1 + } 



(20) 



n=0 p=0 



Therefore the general solution to the differential equation (J7|), for the choice fflUj) . is 
given by 

y = ai?/i(x) + 6ij/ 2 (^) (21) 
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where a\ and b\ are arbitrary constants and y± and y<i are given by (|T9|) and (1201) 
respectively. It is clear that the quantities ?/i and y 2 are linearly independent functions. 
From (IB^|) - (15dj ) and fT2"Tl) the general solution to the Einstein-Maxwell system can be 
written as 

e 2A = 22a 

1 + ax y ' 

e 2v = A 2 y 2 (22b) 

p (b — a)(3 + bx) abx 



C ~ {l + bx) 2 2(1 + bx) 



2 



(22c) 



p (l + as) y (a -6) qfa 

C (1 + bx) y (1 + bx) 2(1 + 6x) 2 1 J 



C (1 + bxf 



(22e) 



The result in fl22al) - fl22ep is a new solution to the Einstein-Maxwell field equations. 
Note that if we set a = 0, f)22ap - fl22ep reduce to models for uncharged stars which may 
contain new solutions to the Einstein field equations (Ha|) - (Hc|) . 



5 Elementary Functions 

The general solution (1211) can be expressed in terms of polynomial and algebraic func- 
tions. This is possible in general because the series f[T^]) and (|2"U1) terminate for restricted 
values of the parameters a, b and a so that elementary functions are possible. Conse- 
quently we obtain two sets of general solutions in terms of elementary functions, by 
determining the specific restriction on the quantity a — b + a for a terminating series. 
The elementary functions found using this method, can be written as polynomials and 
polynomials with algebraic functions. We provide the details of the process in the 
Appendix; here we present a summary of the results. In terms of the original variable 
x, the first category of solution can be written as 



V = 

di{\ + ax) 
+d 2 (l + bx)% 



n+l 



l-(n + l)J2 



4a V (2i-l)(n + i)\ 
^\b^a) (2i)\{n - i + 1)\ 



(1 + bxf 



1 + 



E 



4a 



(2i + 2)(n + z + 
{ n + l )^\b-aJ (n - i)\(2i + 3)\ 



1 + bx) { 



(23) 
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for a — b + a = a(2n + 3) (2n + 1) , where d\ and di are arbitrary constants. The second 
category of solutions can be written as 



y = 



d 3 (l + ax)^(l + bx)% 1 + 



3 ^/ 4a y 



(2i + 2)(n + i)! 



(1 + bxf 



n(n - 1) ^ \ b - a) 



(2i + 3)!(n-i-2)! 



i=l V 7 



(2i- l)(n + i- 2)! 
(2i)!(n-i)! 



(1 + 6x) J 



(24) 



for a — 6 + a = Aan(n — 1), where d% and c?4 are arbitrary constants. It is remarkable 
to observe that the solutions (123!) and (p4l) are expressed completely in terms of ele- 
mentary functions only. This does not happen often considering the nonlinearity of the 
gravitational interaction in the presence of charge. We have given our solutions in a 
simple form: this has the advantage of facilitating the analysis of the physical features 
of the stellar models. Observe that our approach has combined both the charged and 
uncharged cases for a relativistic star: when a = we obtain the solutions for the 
uncharged case directly. 

6 Known Solutions 

It is interesting to observe that we can regain a number of physically reasonable models 
from the general class of solutions found in this paper. These individual models can be 
generated from the general series solution (12"T|) or the simplified elementary functions 
fl23l and fl24l . We explicitly generate the following models. 

6.1 Case I: Hansraj and Maharaj charged stars 

For this case we set a = 0, b = 1 and < a < 1. Then from ( fT9i) we find that 




(y/(l -^)(1 + X)) 6 r- wT~ T [" fj, wT— 

— ~. h... + Vl 1 -a)(l + x) Vl 1 - a )\ l + x ) 



(V(l-q)(l + x)) 3 , (^/(l-aXl + x)) 5 



3! 5! 




= c cos a/(1 — a)(l + x) + c a/(1 — a)(l + x) sin a/(1 — a)(l + x). 
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Equation ([20]) gives the result 



y 2 = co(Vl +x) 3 
3c 



n=0 p=0 



■(1-a) 



(2p + 5)(2p + 2) 



2n+2 



V(l-a)(l + a;) 



3! 



5! 



v/(l-a)(l + x) 



(V(l-«)(l + x)) 2 , (V(l-a)(l + a;)) 



2! 



+ 



4! 



3c 



sin 



(v / r^) 3 l 

Hence the general solution becomes 



- a)(l + x) - a/(1 - q)(1 + s) cos -«)(! + a;) 



+ 



£>i-£> 2 V(l-a)(l + z) cos^/Cl-^Cl + a;) 
D 2 + £>iV(l -a)(l + x)l sin v/(l-a)(l + x) 



(25) 



where .Di and -D2 are new arbitrary constants. The class of charged solutions (I25p is 
the first category found by Hansraj and Maharaj [19]. 

When a = 0, b = 1 and a = 1 we easily obtain the result 



y — a x + &x(l + x) 



(26) 



from (EH). This is the second category of the Hansraj-Maharaj charged solutions. 
We now set a = 0, b = 1 and a > 1. Then from (fl9|) we obtain 



2/1 



Co 



(a-1) 



2(p + l)(2p-l) 



(v/rr^; 



2n+2 



n=0 p=0 

(y/(^-l)(l + x)) 2 (v/(a-l)(l+x)r 
2! 4! 



(v/(a-l)(l + x)) 6 
6! 



V^a-lXl + x) 1)(1 + ar) 



| (y/(^-l)(l + x)) 3 + (v/(a-l)(l + x)) 5 + 



3! 



5! 



c cosh \/ (a — 1)(1 + x) — c \/(a — 1)(1 + sinh a/(« — 1)(1 + x). 
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Equation ([20]) gives the result 



y 2 = c (Vl + xf 
-3c 



n=0 p=0 



(2p + 5)(2p + 2) 



(^TTx) 



2n+2 



V(a-l)(l + x) + 



(V(q-l)(l + x))' 
3! 



(y/( a -l)(l + x)) 5 
5! 



v/(a-l)(l + x) [1 



| (v/(a-l)(l + x)) 2 | (V(a-l)(l+x)) 4 | ^ 



2! 



4! 



-3c 



, - (sinh - a)(l+x) - y/(l - a)(l + x) cosh ^(1 - a)(l +~xj) . 

(va - 1) 

Therefore, the general solution becomes 

y = [z> 2 - D lV /(a- l)(l + x)j sinh y/(a - 1)(1 + x) 

+ [A - D 2V /(a- l)(l + x)j cosh v/(a- l)(l + x) (27) 

where Di and -D2 are new arbitrary constants. This is the third category of charged 
solutions found by Hansraj and Maharaj. 

The exact solutions ( 1251) . ( 126|) and ( |27|) were comprehensively studied by Hansraj 
and Maharaj [19], and it was shown that these solutions correspond to a charged 
relativistic sphere which is realistic as all conditions for physically acceptability are 
met. The condition of causality is satisfied and the speed of light is greater than the 
speed of sound. 



6.2 Case II: Maharaj and Komathiraj charged stars 



If b — 1, then ([23]) becomes 
y = d\{l + ax) 



1 t ^( 4a Y (2i-l)(n + i)! 



i=l 



a) (2i)!(n-i + l)! 



+d 2 (l + x) 2 



3 ^/ 4a y (2i + 2)(n + i + l)! (1 | ^ 



( n + 1 )^V 1 - a 7 (n-z*)!(2z' + 3)! 
for a - 1 + a = a(2n + l)(2n + 3). When 6 = 1 then (J23J) gives 



d 3 (l + ax) 2 (1 + x) 2 



1 + 



' 4a x ' 



n n 



' i=l v 



(2z' + 2)(n + z)! 
ay (2z' + 3)!(n-z'-2)! 



(28) 



;i + x) 



+d 4 



N -A/ 4a Y (2i - l)(n + i - 2)\ , 
l-n(n-l)^( ) ± — „J,\ ; (1 

»=i ^ 



(2i)!(n-i) 



(29) 
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for a — 1 + a = 4an(n — 1). The two categories of solutions (I28I) and (1291) correspond 
to the Maharaj and Komathiraj [17] model for a compact sphere in electric fields. The 
Maharaj and Komathiraj charged stars have a simple form in terms of elementary 
functions; they are physically reasonable and contain the Durgapal and Bannerji [16] 
model and other exact models corresponding to neutron stars. 



6.3 Case III: Finch and Skea neutron stars 



When a = 0, we obtain 



y 



d x - D 2 VT 



x 



cos 



Vi + x + d 2 + d 1 VT+ 



X 



sin 



X 



(30) 



from (1251) . Thus, we regain the Finch and Skea [20] model for a neutron star when the 
electromagnetic field is absent. The Finch and Skea neutron star model has been shown 
to satisfy all the physical criteria for an isolated spherically symmetric stellar uncharged 
source. It is for this reason that this model has been used by many researchers to model 
the interior of neutron stars. 



6.4 Case IV: Durgapal and Bannerji neutron stars 

If we take a = and n = then 2a + b = 0, and we get 

1 3 

y = d\(l + ax) 2 " (5 — Aax) + d 2 (l — 2ax) 2 



from (123]) . If we set a = — \ [ie. b = 1), then it is easy to verify that this equation 
becomes 

y = c 1 (2-x)^(5 + 2x) + c 2 (l + a;)t (31) 

where C\ = and c 2 = d 2 are new arbitrary constants. Thus we have regained 

the Durgapal and Bannerji [16] neutron star model. This model satisfies all physical 
criteria for acceptability and has been utilised by many researchers to model uncharged 
neutron stars. 



6.5 Case V: Tikekar Superdense Stars 

If we take a = and n = 2 then 7a + b = 0, and we find 



y = d 3 (l + ax){l — Tax) 2 + d 4 



1 + ^(1 - lax) - ^(1 - lax) 2 
2 8 



from If we set a = — 1 (ie. 6 = 7) and let x = a/1 — x then this equation becomes 

7 



y = cix(l — -x 2 ) 2 + c 2 



l ~2 X + 24 X 



(32) 
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where c\ = c/382 and c 2 = —d^fi are new arbitrary constants. Thus we have regained 
the Tikekar [21] model for superdense neutron star from our general solution. The 
Tikekar superdense model plays an important role in describing highly dense matter, 
cold compact matter and core-envelope models for relativistic stars. The Tikekar rela- 
tivistic star falls into a more general class of models with spheroidal spatial geometry 
found by Maharaj and Leach [12]; this class can be generalised to include the presence 
of an electric field as shown by Komathiraj and Maharaj [22]. 



7 Physical Analysis 

In this section we demonstrate that the exact solutions found in this paper are physi- 
cally reasonable and may be used to model a charged relativistic sphere. We observe 
from f)22al) and (I22b|) that the gravitational potentials e 2u and e 2A are continuous in 



the stellar interior and nonzero at the centre for all values of the parameters a, b and 
a. From f!22cj) . we can express the variable x in terms of the energy density p only as 



x 



^- {C[2(6 - a) - a]p- x ± y/C 2 [2(b -a)- afp- 2 + 8C[4(6 - a) + ajp" 1 - 4} 



Therefore from ( I22dl) . the isotropic pressure p can be express in terms of p only. Thus 



all the forms of the solutions presented in this paper satisfy the barotropic equation of 
state p = p(p) which is a desirable feature. Note that many of the solutions appeared 
in the literature do not satisfy this property. 

To illustrate the graphical behaviour of the matter variables in the stellar interior 
we consider the particular solution (125]) . In this case the line element becomes 

ds 2 = -A 2 (A - D 2 y/(1 - a)(l + r 2 )) cos ^(1 - a)(l + r 2 ) 

+ (d 2 + D 1 V(l -a)(l + r 2 )) sin - a)(l + r 2 )] ' dt 2 

+ (1 + Cr 2 )dr 2 + r 2 {d6 2 + sin 2 Odtf). (33) 

For simplicity we make the choice A = 1, C = 1, D\ = 1, D 2 = 4 in the metric (|33|) . 
We choose a = | for the charged and we consider the interval < r < 1 to generate 
the relevant plots. 

We utilised the software package Mathematica to generate the plots for p, p, E 2 
and ^ respectively. The dotted line corresponds to a — | and E 2 7^ 0; the solid line 
corresponds to a = and E 2 = 0. In Fig.l we have plotted the energy density on 
the interval < r < 1. It can be easily seen that the energy densities in both cases 
are positive and continuous at the centre; it is a monotonically decreasing function 
throughout the interior of the star from centre to the boundary. In Fig. 2 we have 
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plotted the behaviour of the isotropic pressure. The pressure p remains regular in 
the interior and is monotonically decreasing. The role of the electromagnetic field is 
highlighted in figures 1 and 2: the effect of E 2 is to produce smaller values for p and p. 
From figures 1 and 2 we observe that the presence of E does not significantly affect p 
but has a much greater influence on p closer to the centre. We believe that this follows 
directly from our choice (fTUj) for the electric field intensity; other choices of E generate 
different profiles as indicated in [18]. The electric field intensity E 2 is given in Fig. 3 
which is positive, continuous and monotonically increasing. In Fig.4 we have plotted 
Y on the interval < r < 1 for both charged and uncharged cases. We observe that 
*p is always positive and less than unity. This indicates that the speed of the sound 
is less than the speed of the light and causality is maintained. Note that the effect of 
the electromagnetic field is to produce lower values for and the speed of sound is 
decreased when a^O. Hence we have shown that the solution (T251) . for our particular 
chosen parameter values, satisfies the requirements for a physically reasonable charged 
body. 



p 




0.2 0.4 0.6 0.8 1.0 

Figure 1: Energy density in geometrised units. 




0.2 0.4 0.6 0.8 1.0 

Figure 2: Matter pressure in geometrised units. 
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0< Figure°3: Elecilfic field in geometrised units. 

dp 
dp 
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_j 1 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 L r 

0.2 0.4 0.6 0.8 1.0 

Figure 4: The gradient dp/ dp in geometrised units. 

8 Discussion 

In this paper we have found a new class of exact solutions (12~TT) to the Einstein-Maxwell 
system of equations (f5aT) - (l5dl) . We demonstrated that solutions (1231) and (!2l"l) in terms 
of elementary functions can be extracted from the general class. If a = then we 
regain uncharged solutions of the Einstein field equations which may be new. The 
charged solutions of Hansraj and Maharaj [19] and Maharaj and Komathiraj [17] were 
shown to be special cases of our general class; the uncharged neutron star models of 
Finch and Skea [20], Durgapal and Bannerji [16] and Tikekar [21] are regained when 
the electromagnetic field vanishes. The simple form of the exact solutions in terms 
of polynomials and algebraic functions facilitates the analysis of the physical features 
of a charged sphere. The solutions found satisfy a barotropic equation of state. For 
particular parameter values we plotted the behaviour of p, p, E 2 and * and showed 
that they were physically reasonable. A pleasing feature of our plots is that we can 
distinguish between charged and uncharged exact solutions. The presence of charge 
leads to smaller values of p,p and ^ in the figures generated. This indicates that the 
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presence of charge can dramatically affect the behaviour of the matter and gravitational 
variables. 
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Appendix 



In the Appendix we derive the elementary functions (123]) and (T5 
On substituting r = into equation ([17]) we find that 

4ai(i- 1) - (( a -6) +q) 
(a — b)(2i + 2) (2« — 1) 

If we set a — b + a = 4an(n — 1), where n is a fixed integer and we assume that a ^ 0, 
then c n+ i = 0. It is easy to see that the subsequent coefficients c n+ 2, c n+ 3, c n+ 4, ... 
vanish and equation (1M|) has the solution 



-n(n — 1) 



4a 



(2i - l)(n + i - 2) 



b-a] {2i)\{n-i)\ 
Then from equation (f!5l) (when r = 0) and (1351) we generate 



Cq ,1 < i < n. 



(35) 



Y 1 = c 



i-»(»-i)E(AV (a 7^ w+ \T 2)! ^ 

^\b-aj {2i)\{n-i)\ 



i=l 



(36) 



where a — b + a = 4an(n — 1). 

On substituting r = | into (fTTl) . we get 



a(2i + 3)(2i + l)-((a-6) + a) 

C^_l_^ — ; TTTT" ~ 7~ ~ Ci , ^ U. 



(37) 



(a-6)(2i + 5)(2i + 2) 

If we set a — b + a = a(2n + 3)(2n + 1), where n is a fixed integer and we assume that 
a / 0, then c n+ i = 0. Also we see that the subsequent coefficients c n+ 2, c„ + 3, c n+ 4, ... 
vanish and equation (!37|) is solved to give 



4a V 3(2i + 2)(n + i + l)! 
6^ ) (n+ l)(n-i)!(2* + 3)! C °' 



1 < i < n. 



(3? 
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Then from equations (|T5l) (when r = | ) and (1381) we generate 



(n + 1) 



£ 

i=i 



4a V (2z + 2)(n + i + 1)! 
(n-i)!(2i + 3)! 



(39) 



where a — 6 + a = a(2n + 3)(2rz + 1). The elementary functions (136]) and (l39j) comprise 
the first solution of the differential equation ffTTj) for appropriate values of a — b + a. 
We take the second solution of (ITT]) to be of the form 

Y 2 = [aX - (a - b)]* u(X) 

where u(X) is an arbitrary function. On substituting Y<i into (ITT]) we obtain 

4X [aX - (a - 6)] « - 2 [2aX + (a - 6)] u - [2a - b + a] u = (40) 

where dots denote differentiation with respect to X. We write the solution of the 
differential equation (140"]) in the series form 



u = Y,c n X n+r , C ^0. (41) 

n=0 

On substituting (jUJ) into the differential equation (140]) we find 

oo 

2(a - b)c r[-2{r - 1) + llX^ 1 - (2(a - 6)c n+1 (n + 1 + r)[2(n + r) - 1] 



n=0 



n+r 



-c n [4a(ra + r) 2 - (2a - 6 + a)]) X 



Setting the coefficient of X r 1 in ( 142]) to zero we find 

(a - 6)c r[2(r - 1) - 1] = 0. 

which is the indicial equation. Since Cq ^ and a 7^ b we must have r = or r 
Equating the coefficient of X n+r in (142]) to zero we find that 

4a(n + r) 2 - (2a - b + a) 



0. 

(42) 



■Cr. 



(43) 



^ n+1 2(a-6)(n + r + l)[2(n + r)-l] - 
We establish a general structure for the coefficients by considering the leading terms. 
On substituting r = in equation (143]) we obtain 

Aai 2 - (2a - b + a) 
Cl+1 = (a-6)(2* + 2)(2*-l) Ci ' ( ' 

We assume that a — b + a = a(2n + 3)(2n + 1) where n is a fixed integer. Then 
c n +2 = from (144]) . Consequently the remaining coefficients c n+ 3, c n+ 4, c n+ 5, ... vanish 
and equation (1441) has the solution 



-fn + 1) 



4a V (2i- l)(n + i)! 
(2z)!(n-z + l)! 
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Co, 1 < i < n + 1. 



(45) 



Then from the equations (14T1) (when r = 0) and (145|) we find 

n+l 



w = c 



~t\ b - a J (2i)!(n-i + l)! 



Hence we generate the result 
F 2 = c [oX-(o-6)]3 



n+l 



n 



dE 



4a V (2i-l)(n + i)! 



^\b-aj (2i)\(n - i + 1)\ 



X 1 



(46) 



where a — b + a = a(2n + 3)(2n + 1). 



On substituting r = | into equation (|4*3"I) we obtain 

a(2i + 3) 2 - (2a -6 + a) 
~ (a-6)(2z + 5)(2i + 2) Q ' 



(47) 



We assume that a — b + a = 4an(n — 1) where n is a fixed integer. Then c n „i = from 
( 1471) . Consequently the remaining coefficients c n , c n+ i, c n+ 2, ... vanish and ( 1471) can be 
solved to obtain 



C; 



4a 



3(2i + 2)(n + i)\ 



b-aj n(n-l)(2i + 3)!(n-z-2)! 



c , i < n — 2. 



(4* 



Then from the equations (14T1) (when r = |) and (14*81) we have 



1 + 



n-2 

E 



4a 



(2* + 2)(n + *)! 



^ \ & - a / (2i + 3)!(n-i-2)! 
Hence we generate the result 



X 1 



Y 2 = c [aX - (a - b)]* X? 



1 + 



n-2 

E 



4a 



(2z + 2)(n + z)! 



n ( n ~ !) ^ \b-aj (2i + 3)!(n-i-2)! 



(49) 

where a — 6 + a = 4an(n — 1). The functions ( 1461) and lj4"9l generate the second solution 
of the differential equation (jTTj) . 

The solutions found can be written in terms of two classes of elementary functions. 
We have the first category of solutions 



Y = D x [aX - (a-6)]3 

3 



1 - (n 



n+l 

DE 



4a 



(2i - l)(n + i)\ 



+D 2 Xi 



1 + 



E 



^\b-aj (2i)\(n-i+ 1)! 
4a V (2z + 2)(n + z' + l)! 



X 4 



(n + l)^V & - a / (n-i)!(2z + 3)! 



-X 1 



(50) 
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for a — b + a = a(2n + 3)(2n + 1), where D\ and D 2 are arbitrary constants. In terms 
of x the solution (J5Uj) becomes 



y = di(l + ax) 2 
+d 2 (l + 6x) ! 



n+l 



s / 4a V (2i - l)(n + i)\ , , x , 
l-^ + DE^J (»)l(n-i + l)l (1 + fa) 

(n+l)^V6-a^ (n-i)!(2i + 3)! 1 ; 



(51) 



where d\ = D\\/b and 0*2 = D2 are new arbitrary constants. The second category of 
solutions is given by 



Y = Di[aX -(a -6))] 5 X2 



1 + 



3 ^/ 4a 



(2z + 2)(n + i)! 
n (^- 1 )^V fo - a 7 (2z + 3)!(n-z-2)! 



■a 



1 /- 4a A' (2^ - l)(n + i - 2)! • 

1 — n(n— 1)> I " — — -X 

~^\ b - a J (2i)\(n-i)\ 



X 1 



(52) 



for a — b + a = 4an(n — 1), where Di and D 2 are arbitrary constants. In terms of x 
the solution (152]) becomes 



d 3 (l + ax) 2 (1 + 6x)2 



3 ^/ 4a 



n ( n_1 ) 5 V 6-0 ,/ (2i + 3)!(n-i-2)! 



(2* + 2)(n + *)! 



+ 6x)' 



+^4 



. >A / 4a V (2i-l)(n + i-2)' 
l-n(n-l)^ — - 77 ^A- - f (1 + bxY 

i=i ^ ' 



(2i)!(n-i)! 

where ^3 = D\Vb and dj = L>2 are new arbitrary constants. 



(53) 
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